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Abstract 

Helium—xenon cooled microreactors are a vital technological solution for portable 
nuclear reactor power sources. To examine the convective heat transfer behavior of 
helium—xenon gas mixtures in a core environment, numerical simulations are 
conducted on a cylindrical coolant channel and its surrounding solid regions. Validated 
numerical methods are used to determine the effect and mechanisms of power and its 
distribution, inlet temperature and velocity, and outlet pressure on the distribution and 
change trend of the axial Nusselt number. Furthermore, a theoretical framework that 
can describe the effect of power variation on the evolution of the thermal boundary 
layer is employed to formulate an axial distribution correlation for the Nusselt number 
of the coolant channel, under the assumption of a cosine distribution for the axial power. 
Based on the simulation results, the correlation coefficients are determined, and a semi- 
empirical relationship is identified under the corresponding operating conditions. The 
correlation derived in this study is consistent with the simulations, with an average 
relative error of 5.3% under the operating conditions. Finally, to improve the accuracy 
of the predictions near the entrance, a segmented correlation is developed by combining 
the Kays correlation with the aforementioned correlation. The new correlation reduces 


the average relative error to 2.9% and maintains satisfactory accuracy throughout the 


entire axial range of the channel, thereby demonstrating its applicability to turbulent 
heat transfer calculations for heltum—xenon gas mixtures within the core environment. 
These findings provide valuable insights into the convective heat transfer behavior of a 
helium—xenon gas mixture in a core environment. 
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temperature 


friction factor ò 


radial variation of fully developed 
temperature profile 
convective heat transfer coefficient 


constant in Eq. 10 
eigenvalues 


distance required downstream of 
first measuring point to observe a 


decrease in ANu 


distance required upstream of last 
measuring point to observe an 


increase in ANu 


eddy diffusivity (m7/s) 


(éi K)) u ` viscosity (Pa-s) 

thermal conductivity (W/(m: Ku m circular constant 
hydraulic diameter of coolant channel (m) p density (kg/m°) 

mass flow rate (kg/s) ọ simplified Fn 

Nusselt number m simplified an 

10% of average Nu obtained in each run 

pressure of standard condition (MPa) Subscripts 

Prandtl number 

Peclet number avg average 

power of standard condition (W) b bulk 

heat flux (W/m?) cor calculated by correlation 
Reynolds number in inlet 

dimensionless coolant channel radius m momentum 

inlet temperature of standard condition (K) out outlet 
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velocity (m/s) W 


wall 


1. Introduction 

Nuclear microreactors typically generate up to 20 MW of power and offer 
numerous advantages, including flexible siting, long operational endurance, high safety 
level, and high reliability!!!, They can provide power and heat to both industrial 
facilities and remote communities, and are promising for use in deep-space and deep- 
sea exploration”!. To develop microreactors that can generate megawatt-level electrical 
power, scientists typically employ gas-cooled reactors and the Brayton cycle?!!! to 
satisfy the requirements of small size, light weight, and simple layout. Furthermore, an 
appropriate working fluid must be selected to achieve the desired objectives. 
Incorporating a certain amount of xenon in a helium working fluid can significantly 
enhance its density and compressibility'®!, thus addressing the disadvantages of helium 
while satisfying the aforementioned requirements. Consequently, gas-cooled reactors 
with direct cooling and closed heltum-xenon Brayton cycles have been investigated 
extensively for the development of micronuclear power sources. 

To satisfy the objectives mentioned above, the molar concentration of xenon in a 
helium—xenon gas mixture is typically set between 8.6% and 30.0%. Within this range, 
the Prandtl number (Pr) of the mixture is typically between 0.16 and 0.30!7!, Based on 
reported data, the Pr range of the heltum—xenon gas mixture is lower than that of typical 
gases such as air, hydrogen, and helium"! but higher than those of liquid metals such as 
sodium"! and lead—bismuth eutectic!!”!, This unique Pr range endows the helium—xenon 
gas mixture with convective heat transfer characteristics that are distinct from those of 
the aforementioned fluids. Furthermore, in the core flow channel under conditions of 
high speed, high power, and non-uniform power distribution, the gas properties and 
structure of the thermal boundary layer undergo significant changes, which 
consequently affect the convective heat transfer. The preceding discussion clarifies that 


investigations into the convective heat transfer characteristics of a heltum—xenon gas 


mixture in an intricate core environment is one of the primary areas of research 
pertaining to the thermal hydraulics of helium—xenon cooled microreactors. 

Numerous significant findings regarding the heat transfer characteristics of 
helium—xenon gas mixtures have been published. In terms of experiments, Taylor et 
al.!!] examined the heat transfer properties of mixed gases, including helium—xenon 
gas mixtures, in a uniformly heated circular tube and assessed the feasibility of several 
typically used heat transfer correlations. To examine the convective heat transfer 
characteristics of heltum—xenon gas mixtures with various channel geometries, 
Nakoryakov et al.!!7!'6! investigated the flow and heat transfer properties of a mixture 
in uniformly heated circular and triangular channels, as well as in circular-triangular 
transition regions. To achieve the actual geometric structure of a core, Makarov et a! 
conducted experimental and numerical studies on the heat transfer process of a helium- 
xenon gas mixture in a thin-walled quasi-triangular pipe, which is suitable for densely 
packed fuel elements in the core. They analyzed the effects of channel geometry 
modifications and boundary layer development on the Nusselt number (Nu) and 
discussed the significant decrease in wall temperature at the exit of the quasi-triangular 
reign and its underlying causes. In addition, Qin et al.!'*! examined the convective heat 
transfer of heltum—xenon gas mixtures in a uniformly heated vertical circular tube. They 
developed a more accurate Nu correlation and identified an applicable range. 

Meanwhile, in terms of simulation, Vitovsky et all and Lushchik et al. Di 
conducted researches on the flow and heat transfer properties of helium—xenon gas 
mixtures in a uniformly heated small-diameter (5.5 mm) tube. They proposed a method 
using the mass-average recovery temperature as the characteristic flow temperature, 
which extended the applicability range of the Dittus—Boelter correlation. In contrast to 
the approach of Vitovsky et al., Zhou et al. validated the existing turbulent Pr”! and Nu 


correlations!!! and proposed new correlations with higher accuracy based on numerical 


simulations and theoretical derivations. Additionally, researchers have performed 
numerical simulations to investigate the heat transfer characteristics of helium—xenon 
gas mixtures in various channel types, including trilobe channel?" annular channels, 
narrow rectangular channels'”*!, wire-wrapped annular channels'**!, and channels with 
internal vortex generators!*!, Additionally, to simulate the actual operating conditions 
of a nuclear reactor more accurately, Meng et al Dë! developed a 1/12-scale model of a 
helium—xenon cooled space reactor core and investigated the effects of bypass vessel 
cooling, surface radiation, flow channel blockage, and fuel rods on the heat transfer 
characteristics of the core. 

The main objective of most relevant experimental and simulation studies is to 
establish the Nu correlation, which serves as a critical foundation for the development 


27} 28] and system analysis?) Di codes for nuclear reactors. 


of multiphysical coupling! 
However, these studies are based on the fundamental assumption of sufficiently 
developed temperature and flow fields. Nevertheless, in actual reactors, the power 
distribution is nonuniform, and the axial power distribution approximates a cosine 
function"! Di Consequently, the thermal boundary layer inside the channel cannot 
easily reach a fully developed state'**!, thereby affecting the applicability of the Nu 
correlation in the core channels. To investigate the effect of the thermal boundary layer 
evolution on heat transfer, Sparrow et al.°! derived a correlation for the axial Nu 
distribution of channels under uniform heat flux conditions, including the entrance 
region. The accuracy of the correlation was verified by comparing it with experimental 
results. Siegel et al.!°>! expanded Sparrow’s study and derived a relationship between 
the variation in the wall-to-bulk fluid temperature difference and the axial position 
under an arbitrary axial wall heat flux. Based on the studies of Sparrow and Siegel, 


researchers derived and validated relationships for the axial distribution of Nu under 


various conditions, including parallel plates°!, laminar flow”, and a few simple heat 


flux distributions'**!, 

Clearly, investigations into the heat transfer characteristics of helium—xenon gas 
mixtures under nonuniform power distribution conditions are insufficient. 
Consequently, studies that investigate the effects of key operational parameters, such as 
power, flow velocity, and temperature, on the heat transfer characteristics under a 
nonuniform power distribution are few. Moreover, the applicability of the relevant 
conclusions in a nonuniform power distribution environment requires further 
verification. 

Regarding the derivation of the Nu correlation, the correlation for the axial Nu 
distribution is inadequate when the axial power distribution reflects a cosine function 
and the power density at both the inlet and outlet positions is zero. For microreactors, 
the volume, mass, and complexity of the system should be reduced by eliminating the 
axial reflector, thereby resulting in an axial power distribution that conforms to the 
aforementioned conditions. Furthermore, additional studies are warranted to evaluate 
the applicability of Sparrow’s method to the turbulent heat transfer of low Pr gases in a 
core environment. 

In this study, computational fluid dynamics (CFD) methods are employed to 
examine the convective heat transfer properties of a helium—xenon gas mixture in a 
cylindrical coolant channel located within the core of the Small Innovative helium— 
xenon cooled MObile Nuclear power System (SIMONS)*?!, A three-dimensional 
multiregion conjugated heat transfer model is developed for one-third of the coolant 
channel. Based on this model, analyses are performed to investigate the effects of 
various factors on Nu, including the axial power and distribution, inlet temperature and 
velocity, and outlet pressure. Subsequently, a theoretical correlation is derived for the 
axial Nu distribution in the channel, under the assumption of a cosine distribution for 


the axial power. Using the numerical simulation results, the unknown coefficients in 
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the correlation are determined and a semi-empirical relationship is developed for the 
axial Nu distribution. Subsequently, the accuracy of the semi-empirical correlation is 
tested and compared with those of several existing correlations. Finally, a segmented 
correlation with higher accuracy is obtained by combining the semi-empirical 


correlation with another correlation. 


2. Model and method 
2.1 Geometry model 

The SIMONS energy conversion system is based on a closed Brayton cycle, where 
the core serves as the heat source. The core was designed with a solid structure to 
achieve system miniaturization. Figure 1 shows a schematic illustration of the core. 
The solid core comprised staggered coolant channels and fuel rods. The nuclear fuel 
used was uranium carbide, the moderator was graphite, and the cladding was composed 


of a titanium zircontum molybdenum (Mo-TZM) alloy. 


Coolant 


Flow Heating 
direction section 


Fuel Central axis 
Cladding Reflector ———~ nf --- 
Pressure vessel Adiabatic 
section 


Figure 1. Longitudinal cross-sectional view of the core 


In Figure 1, the dashed hexagonal line on the left represents a radially repeatable 
element located in the core. Because this element exhibits radial symmetry, a radial 1/6 
structure was utilized for modeling. The details of the constructed geometric model are 


shown on the right side of Figure 1. This model comprised solid regions for the fuel, 


moderator, and cladding, whereas the coolant channel was represented by a fluid region. 

Under normal operating conditions, the low-temperature helium—xenon gas 
mixture enters the solid core from the channel inlet at the bottom of the core, is driven 
by the compressor, absorbs the heat released by the nuclear fuel, and flows out from the 
outlet at the top of the core. An adiabatic section is incorporated before the inlet to 
achieve a fully developed flow velocity distribution at the inlet of the heating section. 


The geometric parameters of the model are listed in Table 1. 
Table 1. Geometric parameters of the model 


Parameter Value /mm 
Diameter of fuel 15.0 
Inner diameter of coolant channel 8.0 
Thickness of cladding 0.5 
Pitch 15.0 
Length of coolant channel 1000 


2.2 Thermophysical properties 

For helium—xenon cooled microreactors that utilize a closed Brayton cycle, the 
operating pressure of the helium—xenon gas mixture is typically approximately 2 MPa, 
and the operating temperature is between 400 and1300 KP, To balance the thermal 
power, efficiency, miniaturization, and lightweight requirements of the system, a 
helium—xenon gas mixture with a xenon volume fraction of 12% was selected as the 
working medium for the cycle. 

The temperature-dependent properties of the mixture at a pressure of 1.9 MPa are 
shown in Figure 2. Physical properties, such as specific heat, thermal conductivity, and 
viscosity, are introduced into the CFD software using polynomials that are dependent 
only on temperature. Density was calculated using the ideal gas equation of state. The 
physical properties of the solid regions were simplified to constants, as listed in Table 


2. 


Density /kg/m? 
oo 
Thermal conductivity /W/(m:K) 
Vis cousity /Pa:S 
Specific heat /J/(kg-K) 


0 0 5-0E-06 1078 
500 750 1000 1250 1500 1750 2000 


16 0.4 L1E-04 1084 
12 0.3 8.0E-05 
: 1082 
34 8 02 5.5E-05 ` 
1080 
4 0.1 3.0E-05 
250 


Temperature /K 


—— Density ~~~ Specific heat =—— Thermal conductivity —— Viscosity 


Figure 2. Properties of the heltum—xenon gas mixture over the operation 
temperature 


Table 2. Physical properties of materials 


. : Density Thermal conductivity Specific heat 
ERT, Manel ebe /Wi(m-K) DeK) 
Fuel UC 13630.0 25.3 200.0 
Moderator Graphite 1850.0 70.0 1835.0 
Cladding | Mo-TZM 10220.0 118.0 255.0 


2.3 Numerical models 

This study employs the widely utilized CFD software STAR-CCM-+ to conduct 
numerical simulations of flow and heat transfer in nuclear power systems!) [411], To 
balance computational time and resources, the Reynolds-averaged Navier-Stokes 
method was employed to accommodate turbulence fluctuation terms. A turbulence 
model was introduced to solve the Navier-Stokes equations. Currently, k-epsilon and 
shear stress transfer (SST) k-omega models are widely used in high Reynolds number 
(Re) calculations. The SST k-omega model has been utilized in related studies!?!! 1 
[42] and yielded satisfactory results; thus, it was employed in this study. To enhance the 
accuracy and convergence speed, an implicit coupling solver was employed for both 
the fluid and solid regions. 

Pr; is a dimensionless parameter which affects the turbulent thermal conductivity. 


As mentioned previously, the helium—xenon mixture used in this study is a low Pr fluid 


that possesses properties distinct from those of conventional fluids. A study”! showed 
that the Pr; model developed by Kays'*”! can be used to accurately calculate the Pr: of a 


helium—xenon gas mixture, as expressed in Eq. 1. 
-1 


P = Ce | (CPe,)* f ( - ) Eq. 1 
fr Sen e = e — ex T ————— . 
. 2Pr.b t Prep i i CP ep) Prep 


Pr pe 
Pe, = = Eq. 2 


where Pri» is the Pr, of the bulk region (0.85), C is a constant (0.3), and Pe; is the 
turbulent Peclet number calculated using Eq. 2. Under the conditions investigated, the 
temperature did not significantly affect Pr. 

Using Nu as a metric to quantify the degree of convective heat transfer is a widely 
recognized practice. Nu is expressed in Eq. 3 as follows: 


<ih 


Nu = —, Eq. 3 
Hr q 


where / is the convection heat transfer coefficient, which is expressed as shown in Eq. 


4: 


h= Eq. 4 


Here, tp is the bulk temperature, which can be calculated using Eq. 5 as follows: 
K Í CpputdA 


Eq. 5 
mCp a 


b 


3. Simulation results 

A series of numerical simulations are required to achieve a comprehensive 
understanding into the flow and heat transfer characteristics of helium—xenon mixtures 
in the core environment. In the simulation, the channel inlet was modeled as a velocity 
inlet, whereas the outlet was modeled as a pressure outlet. Additionally, the bottom and 
top surfaces of the solid region were regarded as adiabatic, and the axial cross-sections 
of both the fluid and solid regions were assumed to be symmetrical. The fuel region 


was regarded as a volumetric heat source. In addition, owing to the use of a helium— 


xenon gas mixture and metallic cladding, the effects of heat transfer deterioration and 
power shift resulting from fouling deposition were not considered in the calculations 


[44], [45] 


Nineteen data measurement planes were uniformly distributed along the axial 
direction of the flow channel. Each measurement plane was positioned at a distance of 
0.05 m from one another. No data measurement planes were positioned at the flow 
channel inlet or outlet. This is because the wall heat flux at these two positions was 
significantly lower than that at other locations within the channel. Consequently, 
including these positions in the analysis would yield unrepresentative results and 
extreme variations in Nu. Such variations are not conducive to subsequent analyses or 
interpretations. 

To investigate the variables that affect convective heat transfer and quantify their 
respective magnitudes of influence, a sensitivity analysis was performed under five 
distinct boundary conditions. The operating conditions of the study are presented in 


Table 3, where Run-S was considered as the standard operating condition. 
Table 3. Summary of boundary conditions for all operating conditions 


Run Power Inlet Outlet 
No. Distribution Value Temperature Velocity pressure 
s Cosine Q=3289.5 T=955.0 U=121.9 P=1.9 
[W /K /m/s /MPa 
D1 Reflector at inlet 
D2 Reflector at outlet 
Reflector at in & Q T U P 
D3 
outlet 
D4 Uniform 
Ql 0.25*Q 
Ge Cosine oo T U P 
Q4 1.25*Q 
Tl 0.90*T 
x 
Pa Cosine Q Ke U P 
T4 1.10*T 
U1 0.75*U 
U2 Cosine Q T 0.88*U P 
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U4 1.25*U 


PI 0.75*P 
P2 0.88*P 
P3 Cosine Q T U 1.13*P 
P4 1.25*P 

UQI 0.73*Q 0.75*U 

UQ2 0.88*Q 0.88*U 

UQ3 EE 1.12*Q I 1.13*U P 

UQ4 1.32*Q 1.25*U 


3.1 Model validation 

The model was validated by comparing the simulation and experimental results of 
Taylor et al. !!!!, A schematic diagram of the experimental setup is shown in Figure 3; 
the adiabatic and heating sections of the tube were 328.72 and 352.20 mm long, 


respectively, and the inner diameter of the tube was 5.87 mm. 


Heating section 
Adiabatic section KA A ZA 3 EL A / 


Inlet Outlet 


Ce? it os LI 


Constant power with uniform distribution 


Figure 3. Schematic diagram of Taylor’s experimental section 
To verify the simulation results, data from two runs were selected: Run-689H, 
which featured gas Pr and Re ranges closest to those in the current study; and Run- 
715H, which featured gas mole mass and temperature ranges closest to those in the 
current study. The turbulence model, solver, and Pr; model used in the simulation were 
consistent with those described in Section 2. The boundary conditions for the two 
experiments are listed in Table 4, and a comparison of the simulation and experimental 


results is presented in Figure 4. 
Table 4. Parameter settings of Taylor’s experiments 
Run No. M/g/mol Pr Re Tak Pow/Pa Heat flux/W/m? 
689H 83.8 0.25 52350-87373 295.5 481257 96326 
715H 14.5 0.30 19485-34042 303.0 807381 296622 
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Figure 4. Comparison between simulation and experimental results. (a) Run-689H; 
(b) Run-715H. 


As shown in Figure 4, both experiments showed a significant decrease in the wall 
temperatures near the outlet, which is likely due to the axial heat conduction of the tube 
wall) Under two separate operating conditions, the relative error between the 
experimental and simulated data was less than 5%. The simulated Run-689H and Run- 
715H exhibited maximum relative errors of 4.4% and 4.2%, respectively, whereas their 
average relative errors were 2.8% and 1.9%, respectively, except for two points situated 
near the outlet. 

In addition, both the aforementioned model and the core channel model satisfied 


the criteria for grid independence. Figure 5 illustrates the grid partitioning of the core 
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channel model along the axial direction. The grid was refined in the vicinity of the wall 
within the fluid region to ensure that the center of the first layer of the grid cells was 


positioned within the viscous sublayer. 


4mm 


Figure 5. Detailed grid view of the core flow channel model. 


3.2 Flow and temperature field of standard condition 

To analyze the flow and heat transfer characteristics of the helium—xenon gas 
mixture in the core environment, one must examine the velocity and temperature 
distributions in both the flow channel and the surrounding solid region. Figure 6 shows 
the radial velocity distribution at various axial positions within the flow channel. 

Based on Figure 6, the fluid flowed gradually from the inlet to the outlet and the 
gas expanded owing to the heating of the fuel elements, thus resulting in a velocity 
increase. Under normal operating conditions, the maximum local velocity in the flow 
channel reached approximately 199 m/s, which corresponded to approximately 20% of 
the local speed of sound. Additionally, as shown in Figure 6, a significant velocity 
gradient appeared in the near-wall region. This is attributed to the fluid experiencing 
resistance from the adjacent wall, which resulted in a lower velocity compared with that 


at the mainstream flow region. 


202309.00078v1 


chinaXiv 


Velocity (m/s) 
2.0e+02 


1.6e+02 


1.2e+02 


8.0e+01 


4.0e+01 


RRL 


ee (a) z=25L (b) z = 62.5L (c) z= 100L 
Figure 6. Flow velocity distribution at different axial positions. 


Figure 7 presents the temperature distribution in both the flow channel and the 
surrounding solid region. Figures (a)—(c) illustrate the radial temperature distributions 
at different axial positions, whereas Figure (d) depicts the axial temperature 


distributions of all four regions. 


AAA 


9.5e+02 (a) z=25L (b) z= 62.5L (c) z= 100L 
El 
(d) Axial direction 
Figure 7. Temperature distribution in solid and fluid regions 
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Based on Figure 7, the fuel region demonstrates unsatisfactory thermal uniformity, 
primarily because of its relatively low thermal conductivity. By contrast, the moderator 
and cladding regions exhibited better thermal uniformity. Additionally, a significant 
temperature gradient was observed in the near-wall region of the fluid, where the 
cladding temperature significantly exceeded that of the helium—xenon gas mixture. 

The axial distributions of the radial average temperature for each solid and fluid 
region is illustrated in Figure 8. As shown, the temperature peak points of the three 
solid regions were located in the latter half of the flow channel. Furthermore, the 
proximity of a region to the fluid region corresponded to the proximity of its peak 


temperature point to the outlet. This phenomenon arises because of the gradual increase 


in fluid temperature along the flow channel and the cosine distribution of power in the 
axial direction. Additionally, the helium—xenon gas mixture experienced a rapid 
temperature increase in the middle of the flow channel and a slower temperature 


increase at both ends, which is attributable to the cosine distribution of the axial power. 
1600 
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Figure 8. Axial distribution of radial mean temperature in different regions 


3.3 Effects of operating parameters on Nu 
3.3.1 Power distribution 

To render the current study more applicable to practical operating conditions, the 
axial power distribution for four different operating conditions, namely Run-D1 to Run- 
D4, was designed based on the actual configuration of the axial reflector in the reactor 
core and previous heat transfer experiments. Table 3 presents the boundary conditions 


for each operating condition and Figure 9 shows the normalized power distributions. 
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Figure 9. Power distribution under different conditions 


Figure 9 shows that the addition of the reflector significantly increased the 
normalized wall heat flux density near the inlet of Run-D1 and outlet of Run-D2. 
Additionally, although the normalized wall heat flux distributions for Run-D3 and Run- 
S were almost identical, the former demonstrated a significantly higher wall heat flux 
density near the inlet and outlet of the channel than the latter. 

Figure 10 shows the axial distribution of Nu under different operating conditions. 
Table 5 lists the heat transfer characteristic parameters for each operating condition. To 
analyze the effects of different factors on the axial change rate of Nu, ANu was defined 
as 10% of the average Nu obtained in each run. Additionally, two dimensionless 
distances were introduced: y, which represents the distance required downstream of the 
first measuring point to observe a decrease of ANu, and 6, which represents the distance 
required upstream of the last measuring point to observe an increase of ANu. These 


values are also marked in Figure 10. 
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Figure 10. Axial Nu distribution under different axial power distributions. 


Table 5. Characteristic parameters under different axial power distributions 


Nu : 
SUNS Maximum Minimum Average T g Kee 
S 92.68 44.20 71.73 16.76 4.22 
D1 93.74 48.13 71.59 9.57 5.16 
D2 88.92 56.34 72.96 23.78 17.20 aL 
D3 92.04 56.60 72.62 11.97 = 17.34 
D4 89.34 66.88 172.59 5.81 82.96 


As shown in Figure 10, Nu decreased along the flow direction, and the decreasing 
trend was consistent across all operating conditions, except for the case of uniform 
power distribution. In the case of uniform power distribution, Nu decreased 
significantly near the inlet, whereas the rate of decrease at other positions decelerated 
considerably. Similar phenomena were reported by Sparrow! and Vitovsky!!®. 

Table 5 indicates that the addition of a reflector at the inlet increased the axial 
decline rate of Nu near the inlet. By contrast, incorporating a reflector at the outlet can 
reduce the axial decline rate of Nu, thus improving Nu near the outlet and contributing 
positively to heat transfer. Hence, incorporating reflectors near the core outlet can 
potentially provide supplementary thermal safety margins for the core. However, in the 
case of microreactors, one must scrutinize the tradeoff between the benefits of enhanced 
thermal safety and the associated increase in core volume and weight resulting from the 
addition of reflectors. 

3.3.2 Power value 

To investigate the effect of power on the channel heat transfer, four sets of 
operating conditions, namely Run-Q1 to Run-Q4, with different power values were 
designed, and their Nu distributions were compared with those of Run-S. The boundary 
conditions for each operating condition are listed in Table 3. The axial distribution of 
Nu and the heat transfer characteristic parameters for each operating condition are 


presented in Figure 11 and Table 6, respectively. 
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Figure 11. Axial Nu distribution under different power values 


Table 6. Characteristic parameters under different power values 


Nu 


— Run-Q2 


15 100 


125 


EE Maximum Minimum Average i Ge 
S 92.68 44.20 71.73 4.22 
Ql 78.52 39.97 69.96 3.18 
Q2 85.10 44.34 72.87 3.71 z/L 
Q3 89.03 44.96 72.70 4.04 
Q4 94.22 43.60 70.51 4.44 


As shown in Figure 7, power exerted a positive effect on Nu near the inlet. 
However, as the power increased, the effect weakened. Furthermore, the Nu in the axial 
middle section of the channel decreased gradually as the total power increased. This 
can be attributed to the lower fluid temperature at lower power levels, which results in 
a lower viscosity and higher Re under similar flow velocities. Based on the preceding 
discussion, one can infer that under certain circumstances, such as low-power operation 


or transient conditions, the phenomenon in which the cladding temperature at the 


reactor inlet exceeds the expected value must be considered. 


In terms of the rate of change of Nu, the distribution of Nu under certain working 
conditions exhibited a significant nonlinear trend. Hence, only the rate of change of Nu 


near the outlet is presented in Table 6. As shown in the table, the power level did not 


significant affect the rate of decrease in Nu near the outlet. 


3.3.3 Inlet temperature 

Operating conditions Run-T1 to Run-T4 were imposed to investigate the effect of 
the inlet temperature on the axial distribution of Nu. The boundary conditions, axial 
distribution of Nu, and heat transfer characteristic parameters for each case are 


presented in Table 3, Figure 12, and Table 7, respectively. 
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Figure 12. Axial Nu distribution under different inlet temperatures 


Table 7. Characteristic parameters under different inlet temperatures 


Nu : 
aN: Maximum Minimum Average Y ò SS) 
S 92.68 44.20 71.73 16.76 4.22 
T1 104.67 50.75 82.61 20.79 4.20 
T2 97.48 47.52 76.85 20.13 4.24 z/L 
T3 85.90 41.90 67.25 18.68 4.31 
T4 81.01 39.55 63.27 17.99 4.33 


Based on the results presented in Figure 12 and Table 7, as the inlet temperature 
decreased, the overall Nu in the channel increased, and the rate of increase in the 
average Nu was approximately -0.094/K. This phenomenon can be attributed to a 
decrease in the inlet temperature, which resulted to an increase in the fluid density and 
a decrease in the viscosity, thus causing an increase in Re. 

Although reducing the inlet temperature of the core can improve the overall heat 
transfer efficiency, it may decrease the outlet temperature, which consequently affects 


the thermal-electric conversion efficiency of the system. 


As shown in Table 7, the effect of the inlet temperature on the rate of change of 
Nu near the inlet and outlet was insignificant. 
3.3.4 Inlet velocity 

The results of the simulations for Runs-U1 to Run-U4 can be used to investigate 
the effect of the flow velocity on the axial distribution of Nu based on a comparison 
with the results of Run-S. The boundary conditions, axial distribution of Nu, and heat 
transfer characteristic parameters for each case are presented in Table 3, Figure 13, and 


Table 8, respectively. 
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Figure 13. Axial Nu distribution under different inlet velocities 


Table 8. Characteristic parameters under different inlet velocities 


Nu . 
Rumi Maximum Minimum Average Y ? Um 
S 92.68 44.20 71.73 16.76 4.22 
U1 78.65 36.76 56.80 9.47 4.94 
U2 85.89 40.95 64.46 12.85 4.62 z/L 
U3 97.63 47.25 78.62 30.22 3.95 
U4 103.29 49.14 84.89 47.34 3.67 


Based on Figure 13 and Table 8, Nu increased with the flow velocity in the axial 
middle section of the channel at an average increase rate of 0.468 s/m. This is because 
an increase in flow velocity improves the mass flow rate and reduces the average fluid 
temperature, thus resulting in an overall increase in Re. 


Furthermore, as the velocity increased, the positive effect of high velocity on Nu 


near the inlet weakened gradually, thus resulting in a lower axial decrease rate of Nu 
near the inlet. Additionally, the effect of the inlet velocity on the axial change rate of 
Nu near the outlet was insignificant. 

Based on the preceding discussion, one can infer that increasing the inlet flow 
velocity can enhance the overall heat transfer efficiency of the core. However, the rate 
of improvement in the heat transfer efficiency at the inlet may not be as significant as 
that in other regions. Consequently, this increase in the flow velocity can result in higher 
cladding temperatures at the inlet than anticipated. 

3.3.5 Outlet pressure 

The results of Run-P1 to Run-P4 were compared with those of Run S to investigate 
the effect of the outlet pressure on the axial distribution of Nu. During the modification 
of the outlet pressure, the physical properties of the helium—xenon gas mixture were 
adjusted accordingly. The boundary conditions, axial distribution of Nu, and heat 
transfer characteristic parameters for each case are presented in Table 3, Figure 14, and 
Table 9, respectively. 
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Figure 14. Axial Nu distribution under different outlet pressures 


Table 9. Characteristic parameters under different outlet pressures 


Nu : 
SURAN Maximum Minimum Average Y i Se 
S 92.68 44.20 71.73 16.76 4.22 JL 
P1 75.71 35.39 57.50 16.10 4.38 


P2 84.55 39.93 64.81 16.48 4.32 
P3 98.32 48.89 78.51 21.99 4.22 
P4 105.38 53.03 84.96 24.58 4.18 


As shown in Figure 14 and Table 9, Nu increased with the outlet pressure, and 
the average Nu increased at a rate of approximately 0.031/kPa. The increase in Nu can 
be attributed to the increase in fluid density caused by the increase in pressure, which 
results in an overall increase in Re and an improved heat transfer performance. 

Caution must be exercised when increasing the inlet pressure to enhance the heat 
transfer efficiency in the core because it can affect the volume and mass of the 
compressor, thus potentially affecting the flexibility of the micronuclear reactor power 
source. 

As shown in Table 9, the effect of the outlet pressure on the axial change rate of 
Nu near the inlet and outlet was insignificant. 

3.3.6 Inlet velocity under fixed inlet and outlet temperatures 

One approach to determine the operating range of a closed Brayton cycle nuclear 
power system is to vary the core power while maintaining a constant temperature 
increase in the core and then observe changes in parameters such as the system 
thermoelectric conversion efficiency. Using this method, operating conditions Run- 
UO) to Run-UQ4 were established by adjusting the total power to maintain a stable 
outlet temperature based on operating conditions Run-U1 to Run-U4. The simulation 
results of Run-UQ1 to Run-UQ4 were compared with those of Run-S to evaluate the 
effects of varying the total power. The boundary conditions, axial distribution of Nu, 
and heat transfer characteristic parameters for each case are presented in Table 3, 


Figure 15, and Table 10, respectively. 
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Figure 15. Axial Nu distribution under different inlet velocities with fixed 
inlet/outlet temperatures 


Table 10. Characteristic parameters under different outlet pressures with fixed 
inlet/outlet temperatures 


Nu ; 
Ruay; Maximum Minimum Average R S Un 
S 92.68 44.20 71.73 16.76 4.22 
UQI 76.70 38.19 58.60 13.20 4.72 
UQ2 84.63 41.49 65.14 15.41 4.52 z/L 
UQ3 98.67 47.10 78.25 24.55 4.04 
UQ4 106.78 49.19 84.43 27.86 3.87 


As shown in Figure 15, Nu increased with the inlet velocity, and the average Nu 
increased at a rate of approximately 0.431 s/m, which was lower than the result 
presented in Section 3.3.4. 

As shown in Table 10, similar to the conclusion in Section 3.3.4, as the inlet 
velocity increased, the positive effect of the high velocity on Nu near the inlet gradually 
weakened. However, under the operating conditions investigated in this section, this 
weakening trend became less pronounced owing to the change in the power. 
Furthermore, the inlet velocity did not significantly affect the decreasing rate of Nu 


near the outlet. 


4. Correlation for axial Nu distribution 
4.1 Theoretical derivation 

To depict the axial variation in Nu in the channel in the core environment precisely, 
a theoretical framework that can capture the changing behavior of the thermal boundary 
layer is necessary. Sparrow et al. presented a theoretical derivation based on the theory 
of turbulent boundary layers and eddy viscosity. According to their theory, the wall-to- 
bulk temperature difference along the axial direction of a circuit channel under turbulent 


heat transfer conditions of forced convection with a uniform axial power distribution 


can be expressed as shown in Eq. 6°41. 


= AB? z 
> F, exp - al Eq. 6 


n=0 


qL 


where G(r) denotes the dimensionless wall-to-bulk temperature difference in the 
fully developed temperature field; Fn is the coefficient of the series expansion of the 
wall temperature, which is a negative value that depends only on Re, Pr, and G; B2 is 
the eigenvalue and is related only to Re and Pr. The derivation of Eq. 6 is based on the 
following six assumptions: 

I. The fluid properties are constant. 

II. Compared with the radial heat diffusion, the axial diffusion of both molecular 
and turbulent heat is negligible. 

IH. The mean value of the radial velocity is 0. 

IV. Viscous dissipation is negligible. 

V. The turbulence velocity distribution is fully developed throughout the channel. 

VI. The turbulent Pr can be approximated as 1. 

This equation can be used to precisely calculate the distribution of the axial wall- 
to-bulk temperature difference, including that at the thermal entrance region. 


Additionally, G (rf) can be expressed as 
Giré) = -5 D Pa? 
n=0 


Based on Eq. 6, when z = 0 (i.e., at the entrance), the wall-to-bulk temperature 


difference is 0; when z = © Oe, where the temperature is fully developed), the 
temperature difference is a constant. 

For operational conditions with an arbitrary axial power distribution, the power 
distribution can be regarded as a combination of multiple power steps, with each power 
step affecting the downstream temperature distribution'**!, Based on Eq. 6, if no heat 
flux appears on the wall before position z and a uniform heat flux of dq appears after z, 


then the temperature difference downstream Z can be calculated as 


— tylz =dq- zlean) F, el: ZS (z <Z) Eq. 8 


Integrating the equation above yields the following expression for the wall-to-bulk 


temperature difference at Z under an arbitrary axial power distribution®: 


tw — tplz = — di el 10e |- Ze Eq. 9 


Without axial reflectors, the axial power distribution in the core can be 


approximated using a cosine function as follows: 
1 
q(z) = acos [n(z — zl Eq. 10 
After simplification, it can be expressed as 
q(z) = asin (Tz) Eq. 11 
By substituting Eq. 11 into Eq. 9, the equation for the wall-to-bulk temperature 


difference at Z under a cosine power distribution can be obtained as follows: 


2 
oa v i sin (1z)exp Kb 
DE ha Eege Eq. 12 
tw tylz Zy Fa n ABR Ž q 
n=0 exp RET 


2 
Let ap = sea, Therefore, the temperature difference at Z can be expressed as 
e 
mcos (TŽ) — SG ansin (TŽ) 
~ TTN T n 
=i -2y pg EE T Eq. 13 
w # TEE ™” a2 +1? 
The axial distribution of Nu! under the cosine power distribution is represented 


as shown in Eq. 14. 


Nuts —— 
gl L 
oo Eq. 14 
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In Eq. 14, when Ž approaches 1 m from the left side, the following equation is 


satisfied. 


1 

li dE - a,| = - Eq. 1 
SC fz [eo D sin (1Z)exp (anZ) on iad at 

Thus, in the region near the channel outlet, Nu‘! approaches infinity, whereas Nu 


approaches zero (an> 0, Fn< 0). As Z approaches 0 m from the right side, it satisfies the 


following equation: 


1 
Jim fn [cot (mz) — inthe abl — an} =0 Eq. 16 

Hence, near the inlet, Nu! approaches 0, whereas Nu approaches infinity. 

As shown in Eq. 14, Fy and an are critical parameters that affect the Nu. In Sparrow 
and Siegel’s series of studies!*4! B51, the values of Fn and an were obtained by solving 
the eigenvalue system of the Sturm—Liouville type, which is based on the six 
assumptions mentioned above. In this study, the high flow velocity of the helium—xenon 
gas mixture renders viscous dissipation a significant factor that affects the Nu 
distribution under certain operating conditions. Moreover, the gas is heated vigorously 
in the channel, thus causing changes in the physical properties and a pronounced 
acceleration effect. This renders it challenging to satisfy the conditions of constant 
properties and a fully developed turbulent velocity distribution. Therefore, the six 
assumptions mentioned earlier may not be entirely valid, whereas solving the 
eigenvalue equation to determine Fn and an may introduce a significant deviation. 

In addition, coefficients Fn and an in Eq. 14 are related to the characteristic values, 
which result in an infinite number of unknown coefficients that must be determined, 
thereby rendering the equation impractical for engineering applications. To reduce the 
number of unknown coefficients, Fn and an in Eq. 14 were simplified as @ and o. 
respectively, and existing simulation data were used to determine the values of ọ and o 


for various operating conditions. The subsequent analysis indicates that this 


simplification affected the accuracy of the relationship at a reasonable level. Hence, Eq. 


14 can be simplified as follows: 
wu = 2247 Gafet (nz) -———* __|_w} aa 
pw sin (mZ)exp (wZ) 

Using the MATLAB Genetic Algorithm Toolbox, the axial distribution of Nu for 
all operating conditions in Section 3, which presents axial power distributions in the 
form of a cosine function (excluding Run-Q1 and Run-Q2, which exhibited clear 
nonlinear trends), was fitted using Eq. 17. The variations in -ọ and œ with the average 
Re of each operating condition are presented in Figure 16. The average Re of each 
condition was calculated as follows: 

Rein + R€out _ 2m(Hin + Hout) 


Re = 
E 2 TL MinMout 


Eq. 18 
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Figure 16. Relationship between the coefficients and average Re (dashed line 
represents the function graph of the fitted equation, and the average Re is within the 
range of 5.3 x 10* to 1 x 10°) 


When Pr = 0.264 and the average Re was between 5.3 x 10* to 1 x 10°, the 
relationship between ọ and the average Re, as well as that between œ and the average 
Re, were fitted using power functions. The results are expressed as shown in Eqs. 19 


and 20. 
p = —90.72Reang” Eq. 19 
w = 1075.65Regyg Eq. 20 


Here, the Nu at any axial position Z in the flow channel can be calculated conveniently 
using Eq. 17 . In the following section, Z is uniformly substituted with z in the relevant 
correlations to facilitate understanding. 
4.2 Verification and analysis 

As shown in Figure 17, the axial Nu distribution calculated using Eq. 17 was 
compared with the simulation results for all the operating conditions listed in Table 3, 


where the axial power distribution was cosine (excluding Run-Q1 and Run-Q2). 
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Figure 17. Comparison of calculations using Eq. 17 with simulation results 
(relative error within red and black dashed lines is less than 5% and 10%, 
respectively) 


As shown in Figure 17, the most significant difference was indicated in the region 
near the inlet, with a maximum relative error of 75.3%, whereas the average relative 


error was approximately 5.2%. Approximately 94% of the data points showed a relative 


error within 10%, of which and 78% indicated a relative error within 5%. Moreover, in 
the latter half of the channel (z/L = 62.5), almost all the data points showed a relative 
error within 5%. 

Table 11 presents the boundary conditions of the four newly established cases, 
which are distinct from those in Table 3, used to assess the validity of Eq. 17. 


Table 11. Operating conditions to validate Eq. 17 


Run Power Inlet Outlet 
No. Distribution Value Temperature Velocity pressure 
Vi Q 1.05*T 1.10*U 0.88*P 
v2 Corie 1.10*Q 0.90*T 0.80*U P 
V3 0.90*Q T 1.10*U P 
V4 Q 0.95*T 1.10*U P 


Based on the calculations performed using Eq. 18, the average Re of conditions 
Run-V1 to V4 were approximately 5.9 x 10*, 6.9 x 10%, 8.6 x 104, and 9.3 x 104, 
respectively. 

To evaluate the accuracy of Eq. 17, the calculation error and the errors of several 
existing and widely used correlations were compared. Previous studies!!®) 11} Ui 
showed that the Dittus—Boelter!*®!, Churchill!47!, Kays!**!, and Pickett’! methods are 
more accurate for calculating the Nu of a helium—xenon gas mixture compared with 
other correlations. The expressions and applicabilities of the four relationships are listed 
in Table 12. The accuracies of the four methods above and Eq. 17) were compared, as 
shown in Figure 18. Notably, the Nu values obtained from the correlations in Table 12 


were calculated using the local Re, whereas used Reavg was used in Eq. 17 for the 


calculations. 
Table 12. Nut correlations 
Name Range of Re Range of Pr Correlations 
Dittus-Boelter Re >104 0.7<Pr<160 Nu = 0.023Re°®Pr°* 
Nu = 63 q DUT JfPr 
Churchill Re >104 0.001<Pr < 200 EE (1 + Pr®8)5/6 
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Kays Re >10* 0.5<Pr<1.0 Nu = 0.022Re°®Pr° 
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Figure 18. Comparison of calculations using different equations based on 
simulation results of Run-V1 to Run-V4 (relative errors within red and black dashed 
lines are less than 5% and 10%, respectively). 


As shown in Figure 18, Eq. 17 demonstrated high accuracy in most axial positions, 
with a lower relative error compared with the Dittus—Boelter, Churchill, Kays, and 
Pickett correlations. Most of the observations indicated a relative error of less than 10%, 
except for two points near the inlet. In the latter half of the duct, almost all data points 
indicated a relative deviation of less than 5%. Table 13 lists the maximum and average 
relative errors of each correlation for all the data points under the four operating 


conditions listed in Table 11. 
Table 13. Relative error of calculation results using different correlations 


Dittus-Boelter Churchill Kays Pickett Eq. 17 
Relative error Maximin 140.2% 47.7% 16.1% 554% 60.6% 


Average 55.4% 11.1% 15.4% 103% 5.3% 
As presented in Table 13, the Dittus—Boelter correlation showed the highest 


maximum and average relative errors among the five correlations, thus indicating that 
its predictive performance was the least ideal. However, the accuracy of Eq. 17 was the 
lowest, with only approximately one-half of the error indicated in the Pickett correlation, 
thus rendering it the most accurate among the five correlations. Furthermore, based on 
the calculation results for all four operating conditions, the Kays correlation yielded the 
most accurate prediction for the entrance area (z/L < 18.75), with an average relative 
error of 4.5%. In the outlet area (z/L = 106.25), Eq. 17 yielded the highest prediction 
accuracy, with an average relative error of 2.2%. The results above validated the 
theoretical derivation and confirmed the feasibility of Eq. 17 for turbulent heat transfer 
calculations in a helium—xenon gas mixture in a core environment. 

Furthermore, Eq. 17 showed high accuracies for calculations near the outlet but 
exhibited low accuracy near the inlet, in contrast to other correlations. The improved 
accuracy of Eq. 17, particularly near the outlet, can be attributed to the inclusion of the 
thermal boundary layer changes caused by power fluctuations, which were not 
considered in the other correlations. As the power decreases rapidly near the outlet, 
using a correlation that disregards changes in the thermal boundary layer may result in 
significant distortions in the calculations. Furthermore, consistent with the observed 
low accuracy of the calculations near the inlet, previous studies by Sparrow et al Di 
showed similar challenges in achieving accurate calculations in this region. They 
attributed the cause partially to changes in the fluid properties and the uncertainty of 
eddy diffusivities. However, in the present study, in addition to the reasons above, the 
simplification of formulas in the derivation process and the viscous dissipation caused 
by the high-speed fluid may result in inaccurate predictions near the inlet. 

The preceding discussion highlights that although Eq. 17 provides satisfactory 


accuracy for most axial regions in the channel, further improvements can be realized. 


In this regard, the Kays correlation can be combined with Eq. 17 to derive a new 
segmented relationship, as shown in Eq. 21. To evaluate the effectiveness of the 
approach, Eq. 21 was applied to calculate the four operating conditions specified in 


Table 11, and the resulting calculation accuracy is presented in Figure 19. 
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Figure 19. Comparison of Nu between simulation data and Eq. 21 (relative errors 
within red and black dashed lines are less than 5% and 10%, respectively) 


Eq. 21 employs the Kays correlation to calculate Nu in the vicinity of the inlet (z/L 
< 18.75). In this correlation, Re can be approximated using the principle of energy 
conservation. The Nu values in the remaining regions of the channel were calculated 
using Eq. 17. The data presented in Figure 19 show a significant improvement in the 
prediction accuracy achieved via the optimized correlation (Eq. 21). Specifically, the 
maximum and average relative errors reduced to 13.3% and 2.9%, respectively, while 


a satisfactory level of accuracy was maintained throughout the channel axial range. 


5. Conclusion 

Numerical and theoretical investigations were conducted in this study to examine 
the heat transfer characteristics of a heltum—xenon gas mixture in a cylindrical channel 
within a reactor-core environment. A multiregion conjugated heat transfer model of the 
channel was established, and several special operating conditions were established to 
analyze the effects of these factors on Nu. The correlation for the axial Nu distribution 
in the channel was theoretically derived by assuming a cosine distribution for the axial 
power. Upon confirming the unknown coefficients, the accuracy of this correlation was 
meticulously examined and compared with those of other correlations. The main 
conclusions obtained were as follows: 

(1) The addition of a reflector at the inlet accelerated the axial decline rate of Nu 
near the inlet, whereas the addition of a reflector at the outlet presented the opposite 
effect, which enhanced Nu near the outlet. In the case of the cosine power 
distribution, the Nu near the inlet decreased significantly as the total power 
decreased, whereas the opposite trend was observed in the axial middle section of 
the channel. 

(2) Heat transfer was affected by the inlet temperature, inlet velocity, and outlet 
pressure through their effects on the overall Re of the channel. An increase in the 
inlet temperature resulted in a decrease in the overall axial Nu, whereas an increase 
in the inlet flow velocity or outlet pressure resulted in the opposite trend. Increasing 
the inlet velocity reduced the axial decline rate of Nu near the inlet, whereas 
simultaneously increasing both the power value and inlet velocity alleviated this 
phenomenon. 

(3) A theoretical correlation was established to calculate the axial distribution of 
Nu in a channel with cosine-distributed axial power. The coefficients used in the 


correlation were determined based on simulation data, which resulted in a semi- 


empirical correlation. 

(4) Within the ranges of Re and Pr simulated in this study, the obtained semi- 
empirical correlation demonstrated a high degree of accuracy, with an average 
relative error of 5.3%. Moreover, by integrating the semi-empirical correlation with 
the Kays correlation, the overall relative error was further reduced to 2.9%, and the 
calculation accuracy was satisfactory across the entire axial range of the channel. 

In future studies, a comprehensive approach encompassing both experiments and 
simulations shall be adopted to extensively investigate the heat transfer characteristics 
of the helium—xenon gas mixture. This approach contributes to the verification and 
expansion of the methodology outlined herein. Ultimately, the newly established 
correlation will serve as a foundation for the development of a subchannel analysis code 
designed specifically for helium—xenon cooled microreactors. 
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